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ABSTRACT 



Individuals with spina bifida myelomeningocele (SBM) exhibit brain abnormalities in cortical thickness, white 
matter integrity, and cerebellar structure. Little is known about deep gray matter macro- and microstructure in 
this population. The current study utilized volumetric and diffusion-weighted MRI techniques to examine gray 
matter volume and microstructure in several subcortical structures: basal ganglia nuclei, thalamus, hippocampus, 
and amygdala. Sixty-six children and adolescents (ages 8-18; M = 12.0, SD = 2.73) with SBM and typically de- 
veloping (TD) controls underwent T-p and diffusion-weighted neuroimaging. Microstructural results indicated 
that hippocampal volume was disproportionately reduced, whereas the putamen volume was enlarged in the 
group with SBM. Microstructural analyses indicated increased mean diffusivity (MD) and fractional anisotropy 
(FA) in the gray matter of most examined structures (i.e., thalamus, caudate, hippocampus), with the putamen 
exhibiting a unique pattern of decreased MD and increased FA. These results provide further support that SBM 
differentially disrupts brain regions whereby some structures are volumetrically normal whereas others are re- 
duced or enlarged. In the hippocampus, volumetric reduction coupled with increased MD may imply reduced cel- 
lular density and aberrant organization. Alternatively, the enlarged volume and significantly reduced MD in the 
putamen suggest increased density. 

© 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Spina bifida myelomeningocele (SBM), identified by the herniated 
protrusion of the spinal cord and meninges from the vertebrae at birth 
(Detrait et al., 2005), is the most common congenital birth defect affect- 
ing the central nervous system compatible with survival. It occurs in ap- 
proximately .3-5 per 1000 live births (Au et al., 2010). Almost all 
individuals with SBM incur structural brain abnormalities, including 
variable anomalies of the midbrain and corpus callosum. Most consis- 
tently observed is the Chiari II malformation of the hindbrain, which oc- 
curs in over 90% of cases with SBM (Juranek and Salman, 2010). 
Resultant obstruction at the level of the fourth ventricle leads to hydro- 
cephalus in the majority (roughly 80%) of people with SBM, with many 
receiving shunt diversion as treatment. 

The high comorbidity of hydrocephalus and shunt treatment in indi- 
viduals with SBM disrupts brain macro- and microstructure (Del Bigio, 
2010), leading to variable neural and cognitive outcomes (Hampton 
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et al., 2011). Although intellectual disability is infrequent in SBM, the 
neurobehavioral profile of SBM consists of strengths in lower-order cog- 
nitive and motor domains, such as association and motor learning. 
Weaknesses are in higher-order, contextually based cognitive process- 
es, including problem solving and on-line information integration 
(Dennis and Barnes, 2010). There are often relative strengths in motor 
learning and adaptation to environmental stimuli, but weaknesses in 
dynamic motor regulation such as impaired movement coordination 
and fine motor control (reviewed in: Dennis and Barnes, 2010). 

Given the characteristic features of the Chiari II malformation and 
motor impairments, initial neuroimaging studies of SBM targeted cere- 
bellar abnormalities, revealing generally reduced cerebellar volumes in 
these individuals relative to typically developing (TD) peers (Dennis 
et al., 2004; Juranek and Salman, 2010). However, despite overall vol- 
ume reductions, individuals with SBM show significant volumetric en- 
largement in anterior cerebellar portions and reductions in posterior- 
inferior regions (Juranek et al., 2008). 

At the level of the cerebrum, individuals with SBM also exhibit ante- 
rior to posterior patterning in the cortex as well as in the corpus 
callosum (Hannay et al., 2009; Hannay et al., 2008; Juranek et al., 
2008; Treble et al., 2013). Cortical thickness is significantly increased 
in frontal regions, but decreased in posterior regions (Juranek et al., 
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2008) . When examining cortical complexity (i.e., gyrification) a similar 
pattern emerges whereby frontal regions generally have significantly 
less complexity (except the middle frontal region) despite significantly 
greater complexity in parietal and temporal cortices relative to TD con- 
trols (Treble et al., 2013). Posterior regions of the corpus callosum may 
also be more affected in SBM; the genu is generally well preserved 
whereas the splenium shows the greatest dysgenesis (Hannay et al., 

2009) . 

White matter (WM) and gray matter (GM) volumetric and micro- 
structural alterations also exist in individuals with SBM. WM integrity 
is particularly reduced in long association fiber tracts connecting poste- 
rior and anterior brain regions (Hasan et al., 2008a; Ou et al., 2011). 
Most studies have reported decreased fractional anisotropy (FA) and in- 
creased mean diffusivity (MD) in SBM relative to controls (Ou et al., 
2011; Hasan et al., 2008a; Williams et al., 2013). Lower FA has also 
been observed in the genu of the corpus callosum, but not in the anteri- 
or commissure (Herweh et al., 2009). In a recent probabilistic 
tractography study of attention pathways in SBM, posterior pathways 
showed greater reduction in WM integrity (decreased FA) than frontal 
pathways (Williams et al., 2013). Tectal beaking occurs with high fre- 
quency in people with SBM and was further associated with WM alter- 
ations in both anterior and posterior attention pathways. Limbic tract 
integrity, as previously examined through fiber tracking and region of 
interest (ROI) techniques, is also reduced in SBM, particularly in the for- 
nix and cingulum (Kumar et al., 2010; Ou et al., 2011; Vachha et al., 
2006). 

Little is known about deep GM structure and integrity in SBM. This 
knowledge is relevant given disruptions in periventricular white mat- 
ter, the fornix, and long association fibers, which have connections 
with subcortical GM structures such as basal ganglia nuclei and the hip- 
pocampus. However, no studies to our knowledge have examined spe- 
cific limbic system structures and few have investigated the thalamus 
and basal ganglia structures. Of the three studies that have investigated 
specific basal ganglia nuclei and periventricular GM regions using quan- 
titative methods, the most consistent finding is reduced GM integrity in 
the caudate nucleus of the basal ganglia (Hasan et al., 2008b; Kumar 
et al., 2010; Kumar et al., 2011). However, discrepancies exist between 
these three studies of GM integrity. While decreased FA and increased 
MD values were observed in the caudate nucleus of children with SBM 
relative to controls in one study (Kumar et al., 2010), two other studies 
reported increased FA in this region relative controls (Hasan et al., 
2008b; Kumar et al., 2011). 

Similarly, some discrepancies exist with regard to published reports 
of putamen integrity. One study reported no significant differences in 
SBM compared to controls in DTI metrics of the putamen (Hasan et al., 
2008b) and another study reported increased FA in this structure 
(Kumar et al., 201 1 ). Kumar et al. (201 1 ) also investigated GM integrity 
of the thalamus, but reported no significant differences between pa- 
tients and controls (Kumar et al., 201 1 ), despite its proximal location 
to periventricular WM, which has abnormal features in SBM (Hasan 
et al., 2008a). To date, no studies have examined whether volumetric 
differences in these regions (i.e., limbic, basal ganglia, thalamus) are ev- 
ident in SBM. The clinical implications of elucidating subcortical GM vol- 
ume and integrity differences are relevant since most GM structures are 
clinically read as "normal" during qualitative radiological review 
(Fletcher et al., 2005). 

To further understand subcortical GM in SBM and clarify current dis- 
crepancies in existing literature, we examined the volumetric and mi- 
crostructural properties of deep GM regions in SBM and a typically 
developing comparison group. 

1.3. Hypotheses 

1 Volumetric measures (corrected for total brain volume) of deep GM 
structures will be systematically reduced in participants with SBM rel- 
ative to TD controls. 



2 Measures of microstructure in deep GM will reflect reduced (indicat- 
ed by increased mean diffusivity and reduced fractional anisotropy) 
integrity of these regions in participants with SBM relative to typically 
developing controls. 

2. Materials and methods 

2.1. Participants 

Participants included children and adolescents (N = 66) between 8 
and 18 years (M = 12.0, SD = 2.73) of age with SBM (n = 48) and typ- 
ically developing controls (TD; n = 18). Participants were drawn from 
the Houston cohort of an international project examining neurobehav- 
ioral outcomes of spina bifida (Fletcher et al., 2005). While participants 
with SBM were recruited from clinics at Texas Children's Hospital and 
Shriners Hospital for Children, TD participants were recruited from ad- 
vertisements and community contacts and had no previous or current 
neurological difficulties. Participants with SBM had hydrocephalus 
(20.8% mild at time of testing) and all but one had been treated with a 
shunt near the time of birth. Most participants with SBM had lumbar 
and sacral spinal lesions (83.3%), Chiari 11 malformation (87.5%), callosal 
hypogenesis or hypoplasia (thinning) (95.8%), less than 5 shunt revi- 
sions (81.2%), and infrequent seizure histoiy (none: 66.7%; past: 8.3%; 
current: 4.2%). 

Participants in both groups had verbal and/or nonverbal intelligence 
quotient (IQJ scores between 70 and 120 per the Stanford-Binet Intelli- 
gence Scale: Fourth Edition (Thorndike et al., 1986). Full-scale IQesti- 
mates ranged between 63 and 112 for the group with SBM and 
between 84 and 120 for the TD group. Some children with lower 1Q 
scores were eliminated from the parent study because of concerns 
about their capacity for following directions on neuropsychological as- 
sessment procedures and are not discussed below. All study procedures, 
including MR1 acquisitions, were conducted between 2005 and 2010 
and were in compliance with IRB approval at The University of Houston 
and The University of Texas Health Science Center at Houston. 

Demographic characteristics of the sample are displayed for each 
group in Table 1. As expected, participants with SBM had significantly 
lower verbal and nonverbal 1Q scores (ps < .001 ) than controls and sig- 
nificantly lower socioeconomic status (SES) per the Four Factor Index of 
Social Status (Hollingshead, 1975). The group with SBM also had a 
greater number of non-right handed individuals, which is common in 
samples of SBM (Fletcher et al., 2005) and other samples involving con- 
genital brain injury (e.g., Ross et al., 1987). 

2.2. Magnetic resonance image acquisition 

All MR1 acquisitions were acquired at the University of Texas Medi- 
cal School in Houston using a research-dedicated Philips 3 T scanner 
with SENSE (Sensitivity Encoding) technology and an 8-channel phased 
array head coil. 

22.1. Tj-weighted imaging 

After a conventional scout sequence, high-resolution ^-weighted 
anatomical images were acquired in the coronal plane using a 3D 
turbo fast echo sequence with the following parameters: voxel dimen- 
sions = .94 x .94, slice thickness =1.5 mm, TR = 6.50-6.70 ms, 
TE = 3.04-3.14 ms, flip angle = 8°, DFOV = 240 mm, matrix = 
256 x 256. 

2.2.2. Diffusion tensor imaging (DTI) 

DTI images were acquired in the axial plane using a single-shot spin- 
echo diffusion sensitized echo-planar imaging sequence. Diffusion 
sensitizing gradients were applied in 21 directions (weighting: b = 
1000 s/mm 2 ) with one reference image (b = 0 s/mm 2 ) and the following 
parameters: voxel dimensions = .94 x .94 mm, slice thickness = 3 mm, 
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Table 1 

Demographic data for the spina bifida myelomeningocele (SBM) and typically developing control (TD) groups. 



V al IdUlC 


SBM 


TD 


y vdiutr 


I dll VvlMr 




(n = 42) 


(n = 18) 






Sex (n [% male]) 


24(57.1) 


11 (61.1) 


.421 




Handedness (n [% right]) 


35 (72.9) 


16 (88.9) 


.051 


TD > SBM 


Age at MR1 (M ]SD]) 


12.34 (2.9) 


11.24(2.2) 


.148 




Ethnicity (n [% Hispanic]) 


28 (58.3) 


14 (77.8) 


.282 




Socioeconomic status (M ]SD])* 


31.36(13.2) 


39.83(11.3) 


.018 


TD > SBM 


Full-scale 1Q(M ]SD]) 


85.23 (12.4) 


104.18 (10.3) 


<.001 


TD > SBM 


Verbal IQ(M[SD]) 


85.54 (16.0) 


99.56 (13.0) 


.001 


TD > SBM 


Nonverbal IQ(M[SD]) 


92.00(14.1) 


107.44(14.1) 


<.001 


TD > SBM 



Note: IQand socioeconomic status estimates obtained from the Stanford-Binet Intelligence Scales and Four Factor Index of Social Status, respectively (Thorndike et al., 1986; Hollingshead, 
1975). 



TR = 6500 ms, TE = 65 ms, flip angle = 90°, DFOV = 240 mm, matrix = 
256 x 256. 

23. Magnetic resonance data analysis 

All scans underwent initial radiological evaluation. The basal ganglia 
structures (i.e., caudate, putamen, globus pallidus) were all read as "nor- 
mal" for both groups and the thalamus was read as normal for all partic- 
ipants in the TD group and for 95.8% of participants in the group with 
SBM (2 showed unilateral abnormality and 1 had massa intermedia 
< 1 cm). 

Prior to data preprocessing, anatomical MRI and DTI data were 
reviewed for image quality by an expert user with extensive neuroim- 
aging knowledge and experience (JJ) who was blind to age, gender, 
and neuropsychological performance outcomes. Morphometric analy- 
ses were conducted in FreeSurfer v4.0.5 {http://www.surfer.nmr.mgh. 
harvard.edu) (Fischl, 2012) using a 64-bit Linux computer. Skull- 
stripping and regional segmentation were executed using a fully auto- 
mated command (e.g., recon-all -all) for each Ti -weighted image, yield- 
ing 3 classes of voxels: gray matter (GM), white matter (WM), and 
cerebrospinal fluid (CSF) (Fischl et al., 2002; Fischl et al., 2004a; Fischl 
et al., 2004b). Subsequently, within FreeSurfer, a fully automated rou- 
tine determined delimiting boundaries of deep GM structures (i.e., cau- 
date, putamen, pallidum, thalamus, hippocampus, amygdala). The 
brains of people with SBM are dysmorphic and boundaries are not al- 
ways reliably indicated by the FreeSurfer output. The output images 
from FreeSurfer's segmentation and classification algorithms were visu- 
ally inspected and extensively manually edited by a single rater (JJ) with 
substantial expertise using FreeSurfer's software as well its predecessor, 
Cardviews (Filipek et al., 1989). An example of finalized segmentation 
after manual editing is shown in Fig. 1. At the time of MRI acquisitions 
(collected between 2005 and 2010), FreeSurfer v4.0.5 (released June 
2008) was the most recent version available for analyses. Although sub- 
sequent revisions of FreeSurfer were released, FreeSurfer v5.0.0 (re- 
leased in 2010) was the first version that retained manual edits (i.e., 
persistent edits) and the older version of FreeSurfer was maintained 
to ensure consistency across study subjects. 

For DTI data, within-participant co-registration of the GM masks 
(obtained from the FreeSurfer analyses) in T^space with DTI space 
(B0 or nodif image) were performed using FSLv4.1.0 (FMRIB's Software 
Library, http:/ lwww.fmriblox.ac.uk/fsl) (Smith et al., 2004). Minor head 
motion and eddy currents were corrected in the DTI series with the 
Eddy Current Correction tool included in FSL's Diffusion Toolbox v2.0. 
Skull-stripping and removal of non-brain tissue were performed using 
BET v2.1. FSL's Linear Image Registration Tool v5.5 (FLIRT) (Jenkinson 
et al., 2002) was used to perform a within-participant inter-modal line- 
ar registration between each participant's ^-weighted series (reference 
image) and corresponding non-diffusion-weighted images (b = 0) 
from the DTI dataset. The resultant transformation matrix (and calculat- 
ed inverse) was used to transform the binary segmentation masks of 



each GM structure (obtained from high-resolution ^-weighted images 
using FreeSurfer) to corresponding diffusion-weighted space (detailed 
in Juranek et al., 2012). The diffusion tensors were reconstructed using 
FSL's DTIFIT tool within the Diffusion Toolbox. This transformation ma- 
trix was then applied to each participant to mask individual maps of MD 
and FA for each GM structure. To be as conservative as possible, each 
participant's binary segmentation mask was eroded using a 
2x2x2 mm 3 kernel to reduce contamination from neighboring voxels 
of CSF and/or WM. Final quantitative data were obtained using com- 
mand line utilities to compute average MD and FA circumscribed by 
each eroded deep GM mask in each participant. 

2.4. Data analyses 

Preliminary analyses included examination of the distributions of 
the primary volumetric and DTI variables within and across groups 
and evaluation of hemispheric differences by group in addition to po- 
tential covariates and moderators (e.g., age, sex). IQ was not considered 
as an appropriate covariate. A covariate that is an attribute or an intrin- 
sic characteristic of the condition cannot be meaningfully adjusted to 



B. 




Fig. 1. Examples of final segmentations following extensive manual edits. A) SBM 
participant's Tl -weighted MRI (left panel) with color-coded segmentation boundaries of 
deep gray matter structures (right panel). B) TD participant's Tl-weighted MRI (left 
panel) with color-coded segmentation boundaries of deep gray matter structures (right 
panel). 
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"control" for the attribute, instead creating groups with artificial mean 
differences reflecting the association of the covariate with the disorder 
(Dennis et al., 2009). In addition, handedness was not included as a co- 
variate given the lack of variance in the TD group (Kutner et al., 2005). 

Across groups, distributions of volumetric data were generally nor- 
mal. Within groups, large deviations from group means (~2 SD) were 
rare, but where they occurred, they were specifically evaluated for im- 
pact. The group with SBM had significantly smaller total brain volume 
(TBV) relative to the TD group, F( 1,64) = 24.40, p<. 001. Therefore vol- 
umetric measures were evaluated using corrected ratios (GM structure 
volume to TBV) for all analyses. 

To examine potential hemispheric differences between groups, sep- 
arate 2 (hemisphere: left, right) x 2 (group: SBM, TD) repeated mea- 
sures ANOVAs were examined for each corrected GM volume and DTI 
value (MD, FA). With the exception of thalamus volumes (p = .031 ), 
none of the group x hemisphere interactions were significant for GM 
volumes or DTI metrics (p > .05). Across groups, there were significant 
hemisphere effects (p < .015) found for volumes of pallidum and puta- 
men, FA in all basal ganglia structures (caudate, pallidum, putamen), 
and MD in the amygdala, pallidum, and putamen. However, since 
none of the group x hemisphere interaction effects were significant, in- 
dicating similar lateralization of volumetric and DTI metric values for 
participants in the SBM and TD groups, hemisphere was not included 
in the main analyses below. There were also no significant 
group x hemisphere x sex interactions, indicating hemispheric similar- 
ity for males and females. Given these findings, whole corrected GM 
structure volumes (average of left and right hemispheres) and DTI 
values were computed for each GM structure and examined in subse- 
quent analyses, with the exception of thalamus volumes. 

Across groups, age was not significantly correlated (r < .190) with 
GM volumes, although when groups were examined separately age 
was associated with putamen volumes in the TD, but not in the group 
with SBM. However, age was significantly associated with microstruc- 
ture values both across and within groups separately. Since the litera- 
ture has shown age effects on GM volume and microstructure values 
across development (e.g., Lebel et al., 2008), it was included as a covar- 
iate in all analyses. 

To address the study hypotheses, two sets of results were examined, 
first comparing deep GM volumes and then DTI values between the 
SBM and TD groups on each structure. Each deep GM structure was ex- 
amined using separate ANCOVAs with group (SBM, TD) entered as the 
independent variable. For thalamus volume, for which there was a sig- 
nificant group x hemisphere interaction, a 2-way repeated measure 
ANCOVA with one repeated factor (hemisphere) and one between- 
subjects factor (group) was utilized. For all models, age was included. 




□ SBM 
■ TD 



Putamen 



Hippocampus 



Fig. 2. Group differences on corrected putamen (p = .005) and hippocampus (p < .001 ) 
volumes between participants with spina bifida myelomeningocele (SBM; n = 48) and 
typically developing controls (TD: n = 18). Values reflect averages between left and 
right hemispheres. Note: bars represent standard error. 



3. Results 

3.1. Macrostructure 

Table 2 presents means and standard deviations for each deep GM 
structure volume (with and without correction for TBV, which differed 
between groups). ANCOVA results (with covariate age) indicated signif- 
icant effects (shown in Fig. 2) of group for volumes of the putamen, 
F(l,63) = 8.55, p = .005, rf partial = -119. and hippocampus, 
F(l,63) = 39.73, p < .001, T7 2 p artial = .387, but not for the pallidum, 
amygdala, or caudate, F(l,63) < 1, p > .335. Participants with SBM had 
smaller hippocampi and larger putamen volumes. Age was only signifi- 
cantly positively associated with putamen volume, F( 1 ,63 ) = 5.05, p = 



.028,7] 2 pa r tial = .074. 

For the thalamus, the effects of hemisphere and hemisphere x age 
were not significant (ps > .593), although the hemisphere x group ef- 
fect was significant, F(l,63) = 4.25, p = .043, ff partia i = -063. For the 
between-group effects, neither age nor group was statistically signifi- 
cant (ps > .274). Follow-up examination of the hemisphere x group in- 
teraction revealed a greater, though nonsignificant (ps > .139), 
discrepancy between groups (SBM > TD) in the left (M[SD]: SBM = 
.0063 [.0011]; TD = .0059 [.0005]) relative to the right (M [SD]; 
SBM = .0061 [.0009]; TD = .0060 [.0005]) thalamus. 



Table 2 

Mean (standard deviations) for raw and corrected (raw to total brain volume) deep gray matter structures for participants with spina bifida myelomeningocele (SBM; n = 48) and 
typically developing controls (TD; n = 18). Values reflect averages between left and right hemispheres. 



Brain region 


Raw volume (mm 3 ) 




Corrected volume (mm 3 ) 




SBM 


TD 


SBM 


TD 


Pairwise 
(corrected) 


Total brain volume 


1,183,298.86(97,723.81) 


1,291,996.14(83,270.50) 






SBM < TD *** 


Hippocampus 


4845.14(1060.48) 


7034.00 (812.60) 


.004 (.0009) 


.005 (.0007) 


SBM < TD *** 


Amygdala 


3699.86 (684.60) 


3961.27 (419.55) 


.003 (.0006) 


.003 (.0003) 




Caudate 


7434.11 (1147.82) 


7912.05 (812.64) 


.006 (.0008) 


.006 (.0006) 




Putamen 


12,302.23 (1442.31) 


12,731.45 (1007.44) 


.0104 (.0010) 


.010 (.0006) 


SBM > TD " 


Pallidum 


3682.62 (594.91) 


3944.68 (586.11) 


.003 (.0005) 


.003 (.0004) 


SBM > TB * + 


Thalamus ++ 


14,910.35 (2430.92) 


15,743.18 (1856.24) 


.013 (.0021) 


.012 (.0012) 





Note: reported values reflect the average volume for left and right hemispheres. 
+ Croup differences for pallidum volume were no longer significant after demographic variables (sex, socioeconomic status) were included as covariates. 
++ Between-group differences on thalamus volumes were examined using a 2 (group) x 2 (hemisphere) within-participants ANCOVA. 
* p<.05. 
" p<.01. 
"" p<.001. 
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Table 3 

Means (standard deviations) of DTI values for the spina bifida myelomeningocele (SBM, n = 49) and typically developing control (TD, n = 18) groups. 



Brain region 


Mean diffusivity 






Fractional anisotropy 






SBM 


TD 


Pairwise 


SBM 


TD 


Pairwise 


Hippocampus 


1.089 (.127) 


0.954 (.056) 


SBM > TD *** 


0.178 (.044) 


0.168 (.048) 




Amygdala 


0.894 (.014) 


0.855 (.139) 




0.209 (.030) 


0.186 (.016) 


SBM > TD ** 


Caudate 


0.789 (.049) 


0.751 (.015) 


SBM > TD ■** 


0.163 (.030) 


0.125 (.016) 


SBM > TD *~ 


Putamen 


0.745 (.015) 


0.756 (.011) 


SBM < TD * 


0.158 (.017) 


0.143 (.012) 


SBM > TD *~ 


Pallidum 


0.765 (.029) 


0.774 (.019) 




0.301 (.033) 


0.244 (.027) 


SBM > TD *" 


Thalamus 


0.799 (.035) 


0.775 (.022) 


SBM > TD " 


0.302 (.020) 


0.275 (.014) 


SBM > TD *" 



Note: MD M (SD) reported with values xlO 3 ; reported values reflect the average MD/FA for left and right hemispheres. 
' p<.05. 
" p<.01. 
"" p<.001. 



3.2. Microstructure 

Results for microstructure (MD and FA) are presented in Table 3. 
Only overall models and effects that met critical level of alpha 
(p < .05) are discussed below. 

3.2. J. Mean diffusivity 

There was a significant effect of group on MD values in the thalamus, 
F(l,63) = 9.26,p = .003, ifpjrtia, = .128, caudate, F( 1,63) = 12.57,p = 
•001, rf partial = -166, putamen, F(l,63) = 5.12, p = .027, Tf parUal = .075, 
and hippocampus, F(l,63) = 19.67, p < .001, T) 2 partia i = .238. MD was 
significantly higher for the group with SBM in the thalamus, caudate 
and hippocampus and significantly lower in the putamen relative to 
the TD group. Age was significantly negatively associated with MD in 
the putamen, F(l,63) = 6.08, p = .016, r) 2 part i a i = .088, and pallidum, 
F(l,63) = 21.19,p < .001,77 2 par tiai = -252. 

3.2.2. Fractional anisotropy 

There was a significant effect of group on FA values in the thalamus, 
F(l,63) = 24.34, p < .001, T] 2 partia i = .279, putamen, F(l,63) = 11.22, 
p = .001, ?f partial = .151, caudate, F(l,63) = 25.34, p < .001, impartial = 
.287, pallidum, F( 1 ,63 ) = 39.46, p < .001 , T) 2 part i a i = .385, and amygdala, 
F(l,63) = 10.92, p = .002, T) 2 partial = .148. FA was significantly higher 
for the group with SBM in all structures relative to the TD group. Age 
was significantly positively associated with FA in the pallidum, 
F(l,63) = 4.57, p = .036, rj 2 partial = .068, and in the thalamus, 
F(l,63) = 15.25,p < .001,77 2 partia , = .195. 

3.3. Sex and socioeconomic status as covariates 

3.3.1. Sex 

Given the potential sex differences in CM volumes and microstruc- 
tural metric values, sex was added to the original model as a covariate 
along with age. With covariates sex and age, volumetric results 
remained similar except for reduced group differences in pallidum vol- 
umes (p = .551 ); groups remained significantly different on hippocam- 
pus (p < .001) and putamen (p = .004) volumes. Group differences on 
MD values did not meet the critical value of alpha for the pallidum 
and amygdala, but remained significant in the putamen (p = .047), cau- 
date (p = .001), thalamus (p = .004), and hippocampus (p < .001). The 
addition of sex as a model covariate also did not change results for FA 
values in any region; all deep GM structures, except the thalamus, had 
increased FA in the group with SBM relative to the TD group (p < .001 ). 

3.3.2. Socioeconomic status 

Results were also not substantially altered by the addition of SES as a 
model covariate. When both age and SES were controlled, group differ- 
ences in GM region volume remained similar, with the exception of 
pallidum volume (p = .717). The hippocampus remained significantly 
reduced and the putamen remained significantly enlarged in the 



group with SBM compared to controls. The additional covariate also 
did not substantially alter any of the findings for MD or FA; significant 
group differences in MD persisted in all regions except the amygdala 
and pallidum and groups remained significantly different on FA values 
in all structures except the hippocampus. 

3.4. Examination of shunt status in the group with SBM. 

Given that all but one participant in the group with SBM had re- 
ceived shunt intervention as treatment for hydrocephalus, the number 
of shunt revisions was examined as a predictor for GM volumes and 
MD and FA. Linear regression was used to predict whether more revi- 
sions resulted in smaller corrected volumes and worse integrity com- 
pared to fewer shunt revisions for all examined GM regions. For 
participants with SBM, regression results indicated that number of 
shunt revisions was not strongly associated with GM volumes 
(p > .264), FA values (p > .148), or MD values (p > .080). 

4. Discussion 

Failure of neural tube closure in utero results in an altered pattern of 
programmed fetal development leading to gross anatomical brain 
dysmorphologies. In SBM, the considerable inter-individual variability 
in measures of both structure and function is documented (Fletcher 
et al., 2005). The present study probed distinct structural properties 
(both macro- and microstructural indices) of deep GM in children and 
adolescents with SBM relative to a TD comparison group. 

Contrary to initial expectations of ubiquitous disproportionate re- 
duction in volume and microstructural integrity of subcortical GM 
structures following SBM, findings ( summarized in Tables 2 and 3 ) sup- 
ported a model of regionally specific GM areas being differentially im- 
pacted by this early gestational insult. Whereas the hippocampus 
showed volumetric reduction, the putamen was significantly enlarged. 
The majority of other GM structures (i.e., thalamus, amygdala, caudate, 
globus pallidus) were volumetrically scaled to TBV in both groups. Sim- 
ilarly, results from the DTI metrics (i.e., MD, FA) suggested that deep GM 
microstructural integrity was not uniformly disrupted in SBM. Signifi- 
cant group differences across structures typically followed a pattern 
where the group with SBM showed increased MD values. Three excep- 
tions were in the putamen, where MD was decreased, and in the amyg- 
dala and pallidum, which had similar values between groups. With 
respect to FA values, all evaluated structures (except the hippocampus) 
demonstrated significant group differences with increased FA values in 
the group with SBM relative to controls. 

These findings are consistent with the overall concept that SBM re- 
sults in a complex pattern of disproportionate enlargement in some re- 
gions (e.g., volume of the putamen and anterior cerebellum and 
prefrontal cortical thickness and complexity), disproportionate reduc- 
tion of volume in other regions ( e.g., the hippocampus and inferior-pos- 
terior cerebellum and parietal cortical thickness), in conjunction with 
typical structure of some brain regions (Juranek and Salman, 2010). 
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The anomalous structural patterning in SBM may reflect complex com- 
pensatory mechanisms of neural plasticity in this population. Addition- 
ally, the mechanical effects of hydrocephalus likely complicate the 
protracted brain development in SBM. 

The hippocampus may be particularly affected in SBM given the dis- 
proportionately reduced volume and increased MD observed relative to 
controls. Volumetric reduction coupled with increased MD could imply 
lower cellular density and/or aberrant cellular structure and organiza- 
tion in SBM. The hippocampus develops early in gestation (with individ- 
ual subfields distinguishable around 15-19 weeks) and is particularly 
sensitive to early gestational factors associated with the occurrence of 
SBM, including malnutrition (Arnold and Trojanowski, 1996; Morgane 
et al., 2002) and mechanical injury due to hydrocephalus (Del Bigio, 
2010). Previous studies in clinical and preclinical populations indicate 
thinning and stretching of the fimbria and fornix (Del Bigio et al., 
2003; Naidich et al., 1 982), which likely disrupts subcortical and collosal 
connectivity of the hippocampus (Del Bigio and Zhang, 1998; Mataro 
et al., 2006). Microstructural abnormalities of cellular deterioration in 
axons and synapses in hydrocephalic animal models further support 
the current interpretation of diminished cell density and structural or- 
ganization of the hippocampus in the group with SBM (Kriebel and 
McAllister, 2000). Thus, the macro- and microstructural integrity of 
the hippocampus in SBM may reflect initial damage resulting from fail- 
ure of neural tube closure coupled with prolonged damage resulting 
from hydrocephalus. 

The putamen reflected a unique profile of increased volume, in- 
creased FA, and decreased MD in the group with SBM. The dispropor- 
tionately increased volume and significantly reduced MD values 
observed in the putamen are not indicative of decreased density, but 
rather increased density. The literature of typical development from 
childhood to late adolescence and adulthood supports negative rela- 
tions between age and MD values in the putamen (e.g., Abe et al., 
2008; Lebel et al., 2008; Mukherjee et al., 2002). However, typical devel- 
opment of subcortical gray matter volumes generally peaks in early ad- 
olescence and decreases into young adulthood (Giedd et al., 1999a; 
Giedd et al., 1999b; Good et al., 2001; Ostby et al., 2009). When taken 
together with the larger volumes observed in the putamen, potential 
explanations for group differences in MD include atypical structural or- 
ganization and aberrant apoptosis, or programmed cell death that typi- 
cally occurs across adolescence. Mukerjee et al. (2002) suggested that 
deviations from isotropy could reflect deviations in structural pattern- 
ing during neuronal and glial arborization, which supports our interpre- 
tation of increased FA and decreased MD in the current sample with 
SBM (Mukherjee et al., 2002). However, lack of research in early brain 
structure of infants and young children with SBM complicates elucida- 
tion of group differences. Knowledge of early brain macro- and micro- 
structure in SBM would increase current understanding regarding 
structural and maturational anomalies in SBM. 

Surprisingly, the current analyses also indicated distinct volumetric 
and microstructural patterning between the putamen and caudate in 
SBM. Structural properties of the caudate and the putamen (or dorsal 
striatum) were expected to be similar given the common embryonic 
derivation of these structures from the transient ventricular eminence 
(e.g., telencephalic origin) early in gestation (Marin and Rubenstein, 
2003). In contrast, the globus pallidum emerges from diencephalic ori- 
gins (Miiller and O'Rahilly, 2006) and migrates to its final telencephalic 
position in the median eminence which together with the lateral emi- 
nence forms the floor of the third ventricle (Miiller and O'Rahilly, 
2006). Thus, structural properties of the caudate and the putamen are 
expected to be very similar since neurons within their boundaries are 
derived from the same origin. In contrast, structural properties of the 
pallidum are not expected to share structural similarities with the 
striatal nuclei since those neurons are derived from a different source 
of progenitor cells. Dissimilarity of macro- and microstructure in the 
caudate and putamen may reflect differences in their subcortical loca- 
tion, particularly to the ventricles. Decreased WM integrity within the 



posterior limb of the internal capsule in SBM may result in increased in- 
terstitial space in the extracellular neuropil (Kriebel et al., 1993; Ou 
et al., 2011). The putamen is surrounded by the internal and external 
capsules, which may anchor the putamen and reduce displacement fol- 
lowing hydrocephalus. However, compromised WM integrity of the in- 
ternal capsule in individuals with SBM may allow for the putamen to 
expand. 

While no group differences were observed for the thalamus or cau- 
date volumes, these periventricular structures exhibited a similar pat- 
tern of increased MD and FA in the group with SBM. Given the 
proximity of these structures to the ventricles, increased MD and FA 
may suggest lower density and cellular degeneration following the me- 
chanical effects of hydrocephalus. Periventricular WM disruption 
( stretching) is the most frequent consequence of hydrocephalus and re- 
sults in atrophy and gliosis of periventricular WM, increased extracellu- 
lar space, and neuronal degeneration (Del Bigio, 2010; Del Bigio et al., 
2003; Kriebel et al., 1993). Neuropathological studies in hydrocephalic 
rats indicate that hydrocephalus severity is positively related to thalam- 
ic apoptosis following disrupted periventricular WM (Del Bigio, 2010; 
Mori et al., 2002). Potentially, our observations of increased MD and 
FA in the thalamus and caudate of children and adolescents with SBM 
may be explained by less densely packed cells resulting from induced 
cell death from hydrocephalus. Our finding of increased FA in the cau- 
date in SBM corroborates preliminary findings from a subset of this co- 
hort (Hasan et al., 2008b), which was interpreted as reflecting reduced 
dendritic branching. 

4.3. Limitations 

Functional outcomes across individuals with SBM vary greatly and 
may be complicated by hydrocephalus and associated shunting. As pre- 
viously mentioned, hydrocephalus disrupts neural structure, regardless 
of etiology (Del Bigio, 2010; Scott et al., 1998; Yeates and Enrile, 2005). 
The majority of children and adolescents with SBM in the current sam- 
ple had hydrocephalus and had histories of shunt surgery. The current 
findings of nonsignificant relations between macro- and microstructur- 
al values of GM and number of shunt revisions should be interpreted 
cautiously. Shunting procedures were not uniform across individuals 
with SBM in the current sample and can disrupt brain structure during 
insertion and, when applicable, revision. Both the location and number 
of revisions could have impacted brain structure. Future research could 
benefit from an examination of how, specifically, shunt insertion (and 
revision) alters neural structure in this population and whether shunt 
location influences these outcomes. Therefore, current results represent 
a complex, and potentially synergistic, interplay between neural tube 
defect (i.e., SBM), hydrocephalus (resulting from the characteristic 
Chiari II malformation), and shunting history. 

The current study examined a fairly broad age range, although statis- 
tical attempts to account for age were made. Although age was not 
strongly associated with regional subcortical GM volumes, it was signif- 
icantly associated with the DTI metric measurements in the majority of 
analyses. However, for most models, group accounted for a greater 
amount of variance than age. Additionally, the relations observed be- 
tween age and GM microstructure (negative for MD and positive for 
FA) support previous literature (Lebel et al., 2008; Lenroot et al., 2007; 
Ostby et al., 2009). Methodological weaknesses include the sample 
size of the comparison group and the use of a less advanced neuroimag- 
ing protocol relative to contemporary standards. A larger comparison 
group of TD individuals would be desirable especially given the broad 
age-range assessed. Future research may also benefit from a more ad- 
vanced imaging protocol, including the utilization of greater angular 
resolution by increasing the number of diffusion directions for the DTI 
acquisition. However, we feel that the current data acquisition and pro- 
cessing were sufficient to address the current analytic questions and 
were state of the art at the time of data collection. Lastly, while we did 
not assess for puberty, subsequent developmental changes could have 
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influenced findings, as subcortical GM development peaks after puberty 
(Lenroot et al., 2007). Future research would benefit from an investiga- 
tion into the timing and hormonal effects of puberty on brain matura- 
tion in patients with SBM, particularly given reports of precocious 
puberty in this population (Meyer and Landau, 1984; Oakeshott and 
Hunt, 2003). 

4.2. Implications and future directions 

The current study examined the structural properties of deep GM in 
order to elucidate potential neural correlates underlying disrupted neu- 
robehavioral outcomes. These findings are especially interesting given 
that these areas were largely read as normal in radiological review. Dis- 
ruptions in higher-order cognition and motor functioning following 
SBM may be accounted for by aberrant development of subcortical 
GM. For instance, frontostriatal circuitry, which comprises the prefron- 
tal, basal ganglia, and thalamic regions, is involved in multiple aspects 
of motor learning and control, including motor initiation, inhibition, 
and planning (Alexander et al., 1986). Additionally, the thalamus and 
certain basal ganglia nuclei, particularly the caudate and putamen, are 
shown to facilitate working memory and attentional control, which 
are impaired in individuals with SBM (Baier et al., 2010; Brown et al., 
2008; Burmeister et al., 2005; Cohen and Frank, 2009; Rose and 
Holmbeck, 2007; Vachha and Adams, 2005). Furthermore, limbic sys- 
tem dysfunction, especially hippocampal, in SBM supports diminished 
memory abilities in this population (Dennis et al., 2007; Dennis et al., 
2010; Hampton et al., 2011; Vachha and Adams, 2005; Vachha et al., 
2006). 

Overall, the present results support our concept of regionally specific 
brain areas being differentially impacted by this early gestational insult. 
However, the specific functional impairments associated with such 
structural abnormalities have yet to be identified. Examining neurobe- 
havioral correlates of specific subcortical structures and the underlying 
behavioral differences observed in SBM would promote our under- 
standing of these relations. Furthermore, while informative, macro- 
structural anomalies may only represent CNS damage associated with 
SBM, and do not explain complex interactions between regional micro- 
structure in relation to global brain circuitry. Future research would sig- 
nificantly benefit from longitudinal investigations of subcortical GM 
development in children with SBM. Such investigations may illuminate 
mechanisms of plasticity and provide a foundation for effective inter- 
vention in this population. 
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